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Abstract 

In this paper, we present a unified method for calculating the limit of early exercise boundary 
at expiry. We price American style of general derivative using a formula expressed as a sum of the 
value of European style of derivative and so called American premium. We use the latter expression 
to calculate an analytic formula for limit of early exercise boundary at expiry. Method applied on 
American style plain vanilla, Asian and lookback options yields identical results with already known 
values. Results for selected American style of derivative strategies are compared with limits calculated 
by the PSOR method. 
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1 Introduction 

The growth of variety in financial derivatives traded on markets has increased the need for more general 
and more accurate valuation of their prices. The breakpoint in valuation methods for the financial 
derivatives is dated to early 70. of the 20 th century. The cornerstone laid by Black, Scholes and Merton 
in Black and Scholes (1973) and Merton (1973) or its modifications occur in majority of all pricing 
techniques. The well known Black-Scholes partial differential equation and the theory behind it are 
considered as an important basis in the financial engineering. However, the theory of valuation has 
undergone many changes since that time. 

The most basic classification of the financial derivatives is according to their expiration time (one of 
the main properties). The European style of derivatives can be exercised only at the expiration time T. 
On the other hand, by buying the American style of derivatives the holder obtains a right to exercise it 
at any moment by the expiration time. The early exercise boundary of financial derivative xl — x* (t) 
splits the t — x (time-underlying) space into continuation region C and stopping region S. The derivative 
is exercised if spot value of the underlying is in the stopping region, i.e. (t, xt) £ 5, the derivative is held 
otherwise. 
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We created a unified approach for calculation of the analytic value of limit of the early exercise 
boundary at expiry. Note that position of the early exercise boundary at the expiry was already calculated 
for many particular cases of derivatives Albanese and Campolieti (2006), Alobaidi and Mallier (2006), 
Bokes and Sevcovic (2011), Chiarella and Ziogas (2005), Dai and Kwok (2006), Detemple (2006), Kwok 
(2008), Sevcovic (2008), Wilmott et al. (1995), Wu et al. (1999), etc. In this paper, we present theorem 
that generalizes these methods into one formula. It can be used to determine the limit of exercise 
boundary for American style of a general derivative that can be written in desired form (JTJ) . 

In the second section of this paper, we summarize a method of valuation of the American style of 
derivatives driven by a Brownian motion (it can be used to calculate the value of a derivative in form 
required by following method). In the third section, we present the main result of the paper. The method 
of calculation of the limit of early exercise boundary is presented in Theorem l3.ll Next section consists of 
several examples where the method is applied on plain vanilla options, American style of their strategies, 
on the Asian and lookback options, shout options and British style of vanilla option. In the last section, 
we compare the position of early exercise boundary (of condor spread) calculated by our theorem and 
values calculated by the classical PSOR method. 

2 Value of the American style of derivatives 

In this paper, we analyze the American style of derivative with price given by the formula 



where V eu is the price of European style of derivative, E t is conditioned expected value according to the 
information at time t, I5 is the indicator function for stopping region S and /& is American style bonus 
function. The value of an American style vanilla option in this form was first introduced by Kim (1990). 

Remark 2.1 Price process of the American style derivative discounted by the numeraire is a supermartin- 
gale according to the risk neutral measure. It is the Snell envelope of the pay-off process discounted by the 
numeraire and fl]) discounted by the numeraire is the Doob-Meyer decomposition of this supermartingale. 
For further details see (Karatzas and Shreve, 1988; Chapter 1). 

As an example, we can assume an underlying driven by a stochastic differential equation 



for i G {1, . . . , n} on their domain Del" The values (i' e I, a' > and dW t = (<fW t \ dW t n ) are 
drift, volatility and differential of standard n-dimensional Brownian motion under the joined risk-neutral 
measure Q, respectively. The covariance matrix of dWt is defined for i,j € {1, . . . ,n} by 



where pij £ [—1,1] is the correlation coefficient and pij = 1. 

Let f2 and Af be the pay-off function and the numeraire, respectively. The value V(t,xt) of an 
American style of derivative on the underlying asset Xt is then given by 




(1) 



dx\ = fi % dt + a l dW l tl 



Covar dWi.dWl 



Pijdt, 



V(t,x t ) 



v(t,x t ) + e(t,x t ), 
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where 



and 



v(t,x t ) = N(t,x t )E? {N{T,x T ))~~ n{T,x T ) 



fd(t,x t ) 



e{t,x t ) = Af(t,x t )E? 

B f Q(t,x t } \ p. f Q(t,x t 



l s {u,x u )f d (u,x u )du 



dt 



E 



dx l 



1 ™ 



( Q(t,zt 



i—1 ij — l 

The bonus function in ([1} at expiry is then given by the expression 

f b (T,x T ) = -Af(T,x T )f d (T,x T ). 



dx i dxi 



(2) 



Values of fd on the set of zero measure where / and AT are not differentiable are defined as the arithmetic 
average of limes superior and limes inferior at each point of this set. 



3 Limit of the early exercise boundary at expiry 

In this section, we determine the main result of the paper - the position of the early exercise boundary 
x^p at expiry T for a general class of financial derivatives. The result is stated for a wide class of integral 
equations for pricing American style of derivatives of the form This problem has been already 

considered by many authors for American style of certain derivatives (e.g. Albanese and Campolicti, 
2006; Alobaidi and Mallier, 2006; Bokes and Sevcovic, 2011; Chiarella and Ziogas, 2005; Dai and Kwok, 
2006; Detemple, 2006; Kwok, 2008; Sevcovic, 2008; Wilmott et al., 1995; Wu et al., 1999). Presented 
method is a unified approach solving the generalized problem of finding the position of the early exercise 
boundary at expiry. 

Let KS C K™ be a subset of Euclidean space K™. In what follows, we shall denote by dA a boundary 



of the set A c D with respect to the topology D, i.e. dA = A n D\A 

Theorem 3.1 Consider an American style of derivative V am on the underlying x £ D C M™ with the 
stopping and continuation regions defined by the open sets 5cD and C C D, respectively. Let X£ — 
dS(t, .) = dC(t, .) fort G [0,T] be a (set of) manifold(s) of the early exercise boundary at time t. Suppose 
that the value of V am is given by the equation 



Vam{t,X t ) = V eu (t,X t ) + E* 



l s (u,x u )f b (u,x u )du 



(3) 



where V eu denotes a price of the corresponding European style of derivative and fb(t,x) is a function 
representing the early exercise bonus. Furthermore, we suppose that 

V am (t,x) > fl(t,x) and V am (t,x) > V eu (t,x) for any t € [0, T], x € B, (4) 

where Q(t, x) is the pay-off function at time t for both American style and European style of derivative, 
i.e. 

V am (T, x) = 0(T, x) = V eu (T, x) for any igB. (5) 
Then the limit of early exercise boundary at expiry is given by 

X?r = 8Z+, (6) 

where Z+ = {x T G B; f b (T, x T ) > 0}. 
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Lemma 3.1 Consider a mutually disjoint decomposition D = AL) dAU B of a topological space B, where 
OA = dB. Moreover, consider a set Z so that A C Z C A, where A = A U dA is the closure of the set A, 
then dA = dZ. 



Proof: [of Lemma 13. 1| 

Let a £ dA. For each e > 0, there exists a neighborhood O e (a) so that a £ AnO £ (a) and b £ BPiO e (a). 
This implies, that a £ dZ, i.e. dA C dZ, because a £ Z, but b £ Z C A = A U dA. 

Since dA = A n B\A and ~Z = A, we have <9Z = Z n W\Z C cM and the proof of Lemma follows. □ 



Proof: [of Theorem I3.1| 

Part 1). First, we show that 

5(T,.)c{i T ei;A(r,i T )>0}. 

We have 



T-t 



T 

l s (u,x u )f b (u,x u )du 



' {V am {t,X t )-V eu {t,X t ))>0, 



T-t 



for any t £ [0, T) by (j3]). In the limit t — > T, we can omit the conditioned expected value operator E t 
and we obtain 

ls(T,x T )f b (T,x T )>0. 

If (T,yx) £ S, then we obtain 

fb(T,y T )>0. 



Now suppose that there exists (T, j/t) € 5 such that /h(T, ?/t) = 0. Notice that in the stopping region 5 

a; 



we have the identity V am {t, x) — Q(t, x) for any (t, i) £5 and, consequently, Mj (V am (t, x) — fl(t, x)) = 0. 



Then we have 

d dV dQ, 

= ^ (V am (T, y T ) - n(T, y T )) = -^(T, y T ) - A(T, y T ) - — (T, y T ) 

= l(y eM (T,y T )-fi(T,y T )). 

In the stopping region, exercising the derivative (American style) gives holder higher pay-off than keeping 
it (European style), i.e. 

V eu (t,y) <n{t,y) for (t,y) £S, 

for t sufficiently close to expiry T. According to ([5]), the value of difference between European style of 
derivative and pay-off is increasing (from negative values to zero at maturity). The derivative of this 
difference is positive in the stopping region. This is a contradiction and the proof of first part follows. 

Part 2). Now, we show that 

{x T £ D; f b (T, xt) > 0} C S(T, .) = S(T, .) U X}. 
The function /& can be determined on the stopping region by the following property 
= ^- t (V am (T,y T )-n(T,y T )) 

dV eu _ _ _ C?f2 _ 

= ~^~^ T ^t) - 1 s{T,y T )fb(T : y T ) ~ —(T,y T ). 
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To span the function on the whole domain B, we omit the function lg(u,x u ) and we have 

f b (T,y T ) = ^ (V eu (T,y T ) - n(T,y T )) . 

Notice, that the function fb nullifies movements from the pay-off function Q. 

On the continuous region C, the holder of a financial derivative does not want to exercise it, because 
keeping this derivative yields better pay-off, i.e. 

V eu (t,y)>n(t,y) for {t,y)eC, 

for t sufficiently close to expiry T. According to ([5]), the value of difference between European style of 
derivative and pay-off is decreasing (from positive values to zero at maturity). The derivative of this 
difference is negative and so is the value of function /& in the continuation region. The function fb have 
positive values only on S(T, .). 

Using Lemma I5TT1 we have (JSJ) and the proof follows. □ 

Remark 3.1 Notice, that according to the second part of the proof of Theorem \3.1l we can determine 
function of the American style bonus function fb at expiry by the formula 

fb(T,y)= lim-(V e „(t, y) - Sl(t,y)) . 

t-yT at 

Remark 3.2 The limit of early exercise boundary analyzed in this paper is the expansion of the zeroth 
order. For several financial derivatives of American style, higher order expansion was already calculated. 
Further details on this expansion can be found in Dewynne et al. (1993), Sevcovic (2001), Wilmott et al. 
(1995) for plain vanilla the call option, in Stamicar et al. (1999), Zhu (2006), Zhu and He (2001) for 
the plain vanilla put option and in Bokes and Sevcovic (2011) for average strike Asian options. 



4 Calculation of the early exercise boundary position at expiry 

In this section, we calculate the limit of early exercise boundary at expiry for several types of American 
style of financial derivatives and their strategies. The underlying of all derivatives presented in this section 
is driven by a geometric Brownian motion. Theorem 13.11 does not have limitation on the distribution of 
underlying and can be used also in other models for underlying assets. The method can be used also e.g. 
on Levy processes - the most simple approach is to apply Remark 13 . 1 1 (either analytically or numerically). 

We use parameters r, q and a to denote risk-free continuous interest rate, continuous dividend rate 
and volatility of underlying, respectively. The underlying asset S is driven by stochastic differential 
equation 

dS t = (r-q)S t dt + aS t dW t , (7) 

where Wt is the Wiener process. 

We present examples of application of presented method. Results for the most basic of them are well 
known (e.g. vanilla options or Asian options with arithmetic and geometric averaging). The position 
of early exercise boundary for many other derivatives is calculated by complicated methods or specified 
only by the argumentation without any mathematical formulation. However, we include the derivation 
for comparison purposes. 

Some of the examples presented arc American types of strategies of vanilla options. Trading these 
derivatives is not very common, we use them only to demonstrate Theorem 13.11 on more complex types 
of derivatives. 
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Figure 2: The American style bonus function for a put option with r > q (left) and r < q (right). 



4.1 Plain vanilla options 

The European style of vanilla call/put option gives its holder right to buy/sell the underlying S at 
maturity time T for the expiration price X. The position of the early exercise boundary at expiry for 
call and put vanilla option is equal to boundary of set of positive values of f£ al1 (see Figure [1} and /f"* 
(see Figure^), i.e. 



dZ, 



-{-call 



r 

X,-X 

q 



= S* T cal1 and dZ+ put = min 



r 

X, -X 

q 



T > 



respectively. This result is well known and can be found also in Albanese and Campolieti (2006), Detemple 
(2006), Kwok (2008), Wilmott et al. (1995) and many other sources. 

We derive the position of the early exercise boundary according to the method presented in this paper. 
The pay-off functions for call and put options are 

n call (t,S;X) = (S~X) + and ft put (t, S; X) = (X - S) + , 

respectively. The value of European style of vanilla option (the well known solution of Black-Scholes 
partial differential equation extended by Merton) for both call and put option is given by 



where cL 



C eu (t, S; X) = e -« (T- * ) S* (d t ) - e- r(T - t) X$ [d t - aVT~tj , 
P e „(i, S- X) = e- r{T -^X$> (-d t + aVT~~t) - e -«( T -*)S$ (-dt) , 

_ \n% + (r-q+4-)(T-t) 



(8) 
(9) 



We know the value of both pay-off function and European style of option, so we can calculate /& at 
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^■condor ^-condor ^condor 




Figure 3: The American style bonus function for condor spread with — A4 + A3 + A 2 — X\ > and 
r(X 3 + X 2 -X 1 ) > <?A 4 (left), -X 4 + X 3 + X 2 - X x > and r {X 3 + A 2 - X x ) < <?A 4 (center) and 
-X 4 + X 3 + X 2 - X x < (right) 



the expiry according to Remark 13. 11 i.e. for call option and put option we have 

(I for S < X, { rX - qS for S < X, 



respectively. 



f (9- 



r) for S = X, and fP ut (T, S) = { f (r - g) for S = X. 



qS - rX for S > X 







for S > X, 



4.2 Option strategies 

There are many option strategies that can be considered. European style of these financial derivatives is 
linear combination of plain vanilla options, therefore the derivation is similar for all of them. The position 
of early exercise boundary at expiry is calculated only for very few of them. The results for American 
style strangle spread was presented by Alobaidi and Mallier (2006), Chiarella and Ziogas (2005). 

The most complex frequently used strategy consisting of vanilla options is European style of condor 
spread and its restriction butterfly spread. The European style of condor spread is a linear combination 
of four vanilla call options. There are three different cases, when determining the position of the early 
exercise boundary at expiry by the boundary of set of positive values of f^ ondor _ 



In the first case, if we have — X4 + A3 + X 2 — X x > and r (A3 
boundary points has only one element (see Figured] (left)). 



A 2 — Ai) > (7A4, then the set of 



dZ+ 



condor 





r 




max 


Ai, ~X X 


,A 2 




q 



_ Q*condor 

iJrp 



In the second case, we have — A4 + A3 + X 2 — X\ > and r (A3 
boundary points has three elements (see Figure [3] (center)). 



Q Z +condor 



A 2 - Ai) < qX±. The set of 







r 


,A 2 


| min 


max 


Ai, —X x 

q 







X3, - (A3 
9 



■ A 2 — Ai ) 



, Xa 



S, 



*condor 



In the last case, we have — A4 
(see Figure[3] (right)). 



A3 + A 2 — Ai < and the set of boundary points has two elements 



Q^+condor 







r 


,A 2 


■|min 


max 


Ai, —X x 

q 







^3? _ (^3 + X2 — X\ 



, Xa 



Q*con 



dor 
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We present derivation for the position of the early exercise boundary at expiry for a condor spread 
according to the method presented in this paper. The pay-off function of a condor spread is 

tl condor (t,S;X 1 ,X 2 ,X 3 ,X i ) = Sl call (t,S;X 1 )-n call {t,S;X 2 ) 

-n cal1 (t, S;X 3 ) + n cal1 (t, S; X A ) 

= (S - Xt)+ - (S - X 2 )+ -(S- X 3 )+ + (S- X 4 ) + , 

for Xi < X 2 < X 3 < X4, where the case X 2 — X 3 is called a butterfly spread. The price of a European 
style of condor is calculated by the formula 

Veu nd ° r (ti S'i Xi, X 2 , X 3 , X 4 ) = C eu (t, S; X\) — C eu (t, S; X 2 ) 

—Ceu{t, S; X 3 ) + C eu {t, S; X4), 

where the function C eu is defined by ©. 

Once more, we use Remark 13. II to calculate the bonus function for American style of condor spread. 



fH ondor (T,S) = <( 





^ (q - r) 

qS — rXi 

^L(q + r)- rX x 

r(X 2 -X 1 ) 

^(r-q) + r(X 2 -X 1 ) 
r(X 3 + X 2 -X 1 )-qS 



for S < X u 
for S = X x , 
for Xx < S < X 2 , 
for S = X 2 , 

for X 2 < S < X 3 or X 2 = S = X 3 , 

for S = X 3 , 

for X 3 < S < X^ 



r(X 3 + X 2 - Xx)- ^(q + r) for S = X 4 , 



t r(-X 4 + X 3 +X 2 -Xx) 



for X 4 < S. 



4.3 Asian and lookback options 

Options with some extra properties are usually called exotic options. Subgroup of exotic options called 
path-dependent options contains also Asian and lookback options. Asian and lookback options depends 
on some type of average (A) and the extreme value (maximum M or minimum m) of the underlying 
through the lifespan of the option, respectively. 

We calculate the value of limit at the expiry for American style of floating strike of both Asian and 
lookback options. The pay-off functions of floating strike Asian call, Asian put, lookback call and lookback 
put options are 



04lf n (*. S, A) = (S- A) + , Q±i an (t, S, A) = (A- S) + 



lookback 
call 



(t, 5, m) = (S — my 



put 

and n l p °° t kback (t, S, M) = (M - 5) H 



respectively. Now, we use ([2]) to determine the bonus function fi,. According to Hansen and Jorgensen 
(2000) we can use the numeraire Aft — e rt . We define the general exponentially weighted average A at 



time t by the formula 



A 



e -Ht-u)( Su y dU) 



1 -e jo 

where A = for regular and A > for weighted average. By setting the parameter p to values 0, 1, —00 
and 00, we can construct expressions for continuous geometric average, continuous arithmetic average, 
minimum value and maximum value, respectively. 
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The function fd for call option and for put option has form 



f c d all (t,S,A) = 



for S < A, 

\ (limsu Ps ^ A fX l % S, A) + liminf s _ A /-»(t, S, A)) for S = A, 



e~ rt (- r (S-A)+ns-HA) 



for A < S 



and 



e- rt (-r(A-S)+^ A -^s) for S < A, 

fT (t, S,A)={ I (lim sup 5 _^ f d ut (t, S, A) + lim inf s ^a fT (*, S A)) for 5 = A, 
for A < S, 

respectively. The value fis = ( r — q)S according to J7]) and ^ is drift of the stochastic differential 
equation 

dA = fJL A dt + a A dW t A . 

According to Hansen and Jorgensen (2000) and Bokes and Sevcovic (2011) we have the form of stochastic 
differential equation for geometric average 

1 A 9 
dA 9 = — AS In -4- dt, 



for the arithmetic average 



dA a t =- t {S t - A a t ) dt, 



and for both maximum and minimum the equation has form 



dA; 00 = dA™ = 0. 



Finally, the stochastic differential equation for the exponentially weighted general average has form 

A 



dA t = At 



The bonus function at the expiry for call and put Asian options with continuous geometric average has 
form 

for S < A, 

4 (g - r) for S = A, 



f° all ' 3 (T,S,A) 



and 



fr t ' 3 (T,S,A)={ 



-r A + qS - ±A In | for A < S 



rA-qS + ± A In 4 for 5 < A, 



T S 

4 (r - q) for S = A, 

for A < S, 



respectively. The boundary of set of positive values is given by 



dz +caii, g = J( S ,A)eWI$.;^= max 



S 



and 



dZ +put,g = ^ A) e 



>■+> 



A 



S 
A 



1,G 



1,G 



-y*put,g 



where G is the positive solution of transcendental equation 

r-gG--lnG = 0. 
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The solution G is unique on R + for q > and T > 0. 

The bonus function at the expiry for call and put Asian options with continuous arithmetic average 
has form 

for S < A, 

4 (q - r) for S = A, 

(q+±)S-(r+±)A forAKS' 



f™ ll ' a (T,S,A) = { 



and 



fl ut ' a {T,S,A) 



(q+±)S- 



4(r-«) 



±)A for S < A, 
for 5 = A, 
for A < S, 



respectively. The boundary of set of positive values is given by 



dZ. 



-\-call,a 



(S,A) eR};- = max 



and 



A 

S_ 

'A 



r + i 



9+^ 



1. 



T 



SYrp 



p^*put,a 



q+ t- 

Finally, the bonus function at the expiry for call and put lookback options has form 



/r™(T,S,m) 



for S < to, 

f (q - r) for S = to, 
— rm + qS for to < S 



and 



fr x (T,S,M) 



rM - qS for 5 < M, 
f(r-q) for 5 = M, 
for M < S, 



respectively. The boundary of set of positive values is given by 



and 



dZ+ mtn 



dZ +max 



-{ 



(S,m) € Ri; — = max 
T to 



(5,M)eR^;|: = 



1,^ 

q 



The same results for the limit of early exercise boundary at expiry can be found in Detcmple (2006), 
Wu et al. (1999) for geometric average Asian options, in Sevcovic (2008) for arithmetic average Asian 
options, in Dai and Kwok (2006) for arithmetic average Asian and lookback options and in Bokes and 
Sevcovic (2011) for both geometric and arithmetic Asian options. 

This is the first time that the position of early exercise boundary at expiry for general average is 
mentioned in the literature. It can be derived similarly as the previous examples. The bonus function at 
the expiry for call and put Asian options with general exponentially weighted average has form 



f£ alLw (T,S,A) = < 



for S < A, 

i(q-r) for S = A, 

-rA + qS + A^p^^y-l) for A < S 
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and 



f^(T,S,A) = i A( r -q) 



rA -qS~ A-p^ 

^ p(l — e 



-r^ \ \a) 



respectively. The boundary of set of positive values is given by 



dZ, 



-\-call,w 



(S,A)e 



and 



5 
A 



A 



1,Y 



where Y is the positive solution of transcendental equation 

A 



1 ) for S < A, 
for S = A, 
for A < S, 



r — qY — 



p(l - e- XT ) 

The solution F is unique on R + for q > and T > 0. 



(Y p - 1) = 0. 



4.4 Shout options 

Shout options are financial derivatives similar to European plain vanilla options. The difference is that 
the holder of a shout option can once during the life of derivative "shout" to the writer, i.e. the option 
expires and the strike price is reset to actual spot price of the underlying asset. The shouting action is 
conditioned by in-the-money position of the option. According to this property, we need to know optimal 
shouting boundary along with the limit of the boundary at the expiry. The position of optimal shouting 
boundary at expiry, i.e. the boundary of set of positive values for call shout option is the same as for put 
shout option: 



-shout 



X = S. 



* shout 



This result can be also found in Alobaidi et al. (2011). However, the value of this limit was derived only 
by argumentation and without any mathematical formulation. 

Now, we present the derivation according to the Theorem 13.11 The pay-off function of call and put 
shout option is 

f forS'^A, 
[ S - X + C eu {t,S;S) forS>A 



9. 



call, shout 



(t,S;X) 



and 



shout 



(t,S;X) = 



X - S + P eu (t,S;S) for S < X, 
for S > A, 

respectively. The functions C eu and P eu are defined by © and (O, respectively. 

Notice, that the underlying S is under the same measure as for the vanilla option, thus we have the 
numeraire Aft = e rt . In this case, we use the idea from © to determine the bonus function 



/call, shout frp q\ 
b y 1 ^) 



for S < A, 
oo for S > X 



and 



[■put, Shout I rp Q\ 



oo for S < A, 
for S > X. 
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4.5 British vanilla options 

The British vanilla option is financial derivative hedging the real trend of the underlying asset. This 
feature allows its holder to exercise the option prior to the expiry T and receive the best prediction of 
the pay-off according to the real trend of underlying restricted to the contract drift fj, c . 

The position of the early exercise boundary for British plain vanilla call and put option is equal to 
the boundary of set of positive values of f^ B ' cal1 anc j f^ B ' put , i.e. 



^ry+GB.call 

oZ T = max 



r 

X, X 

q + 



_ Q*GB,call 

tjrp 



and 

oZ T = mm A, A = b T 



■ ^ g*GB ,call 



q + He 

respectively. This result is consistent with the calculation presented in Peskir and Samee (2008a, b). 
However, the derivation presented in this paper is much shorter and straightforward. 

Now, we present the derivation according to the method presented in this paper. The pay-off functions 
of call and put British vanilla option are 

^GB.call = e MT-t) S Q (^e) _ X q ( d »c _ a ^T~~t} 

and 

Q GB,put = x ^ J_ d Mc + a ^/Y~tj - e^-^S® (-«#") , 
respectively. The function <£>( • ) is standard normal cumulative distribution function and the expression 

jMc _ In f + (Mc+^ 2 )(T-t) 

Notice again, that the underlying S is under the same measure as for the vanilla option, thus we have 
the numeraire Aft = e rt . We can use the idea from © to determine the bonus function 

rGB.callfrr m _ / for S < X, 

h ^ S >-^ . ■ - rX ioiS>X 




and 



f GB, P ut (T s)= I rX - (l + Hc)S for S < X, 
b I for5>X. 



5 Verification by numerical method 

In this section, we compare and verify analytic values of the limit of early exercise boundary with values 
calculated by the PSOR (projected successive over relaxation) method introduced in Elliot and Ockcndom 
(1982). For further details on the PSOR method see Kwok (2008) or Wilmott et al. (1995). Values for 
plain vanilla options, Asian options and lookback options were derived in other sources. We present only 
the results of comparison for condor spread as a case with several different values of the limit at the 
expiry (results for all of the other examples are very similar). 

In Table [TJ we present values of the limit of early exercise boundary at expiry calculated by the PSOR 
method with 250 time and 40000 space steps. The space domain was bounded to the interval [—1.5, 1.5]. 
The SOR was calculated with parameter to — 1.4 and tolerance e — 10~ 14 . The expiration time has to 
be chosen close to zero, we set the value T = 10 -8 in the calculation. 

Each numerical value tends to analytic value calculated by the formulae presented in section 2J The 
numerical results are improving with increasing density of the time-space mesh and decreasing tolerance 
(as expected). The relative error of results with the highest precision is of the order 10~ 4 . 
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Table 1: Comparison of analytic and numerical values of the limit of early exercise boundary for condor 
spread. 
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V 

Al 


A 2 


A 3 


A 4 


a theor 




error 


3% 


2% 


30% 


1 


3 


4 


5 


1.5 


1.50002 


0.001% 
















1 


1.00033 


0.033% 


2% 


3% 


30% 


1 


3 


4 


5 


4 


4.00074 


0.019% 
















5 


5.00159 


0.032% 


3% 


2% 


30% 










1.5 


1.50006 


0.004% 


1 


2 


3 


4.5 




















4.5 


4.49844 


-0.035% 



6 Conclusions 

In this paper, we presented a unified method for calculation of a limit of the early exercise boundary 
at expiry. The method is applicable for American style of a wide range of financial derivatives. We 
calculated and compared the analytic value of early exercise boundary limit for plain vanilla options 
as well as for Asian and lookback options. The results coincides with the well known values that can 
be found in Albanese and Campolieti (2006), Alobaidi and Mallier (2006), Bokcs and Scvcovic (2011), 
Chiarella and Ziogas (2005), Dai and Kwok (2006), Detemplc (2006), Kwok (2008), Scvcovic (2008), 
Wilmott et al. (1995), Wu et al. (1999) and in many other sources. Moreover, we also calculated the 
analytic value of the limit for shout option and American style of exponentially weighted floating strike 
Asian option with general average, British vanilla option and option strategies (represented by condor 
spread). We verified values by their comparison with values calculated by the PSOR method. Presented 
straightforward method is simple and unifies the approach to methods so far used for the calculation of 
limit of exercise boundary. 
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